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ABSTRACT 

1 We study the propagation of the ionizing radiation emitted by the secondary 

star in 77 Carinae. We find that a large fraction of this radiation is absorbed by the 
primary stellar wind, mainly after it encounters the secondary wind and passes 
through a shock wave. The amount of absorption depends on the compression 
<^ ■ factor of the primary wind in the shock wave. We build a model where the 

compression factor is limited by the magnetic pressure in the primary wind. We 
find that the variation of the absorption by the post-shock primary wind with 



CN ■ orbital phase changes the ionization structure of the circumbinary gas and can 

account for the radio light curve of rj Car. Fast variations near periastron passage 
are attributed to the onset of the accretion phase. 

cn 

Subject headings: (stars:) binaries: general— stars: mass loss— stars: winds, 
outflows— stars: individual (rj Carinae) 

a: 

P : 1. INTRODUCTION 

The 5.54 yr light periodicity of the massive binary system r] Car is observed from the 
radio (Duncan & White 2003; White et al. 2005) to the X-ray band (Corcoran 2005). 
According to most models, the periodicity follows the 5.54 years periodic change in the 
orbital separation in this highly eccentric, e ~ 0.9, binary system (e.g., Hillier et al. 2006). 
The X-ray cycle presents a deep minimum lasting ~ 70 day and occurring more or less 
simultaneously with the spectroscopic event (e.g., Damineli et al. 2000), defined by the fading, 
or even disappearance, of high-ionization emission lines (e.g., Damineli 1996; Damineli et al. 
1998, 2000 Zanella et al. 1984). The spectroscopic event includes changes in the continuum 
and lines (e.g., Martin et al. 2006, a,b; Davidson et al. 2005). The X-ray minimum and 
spectroscopic event are assumed to occur near periastron passages. 

The radio emission comes from thermal bremsstrahlung (free free) emission from ionized 
material (Cox et al. 1995; Duncan & White 2003). We follow previous suggestions that most 
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of the UV radiation that ionizes material at large distances from the binary system originates 
in the hot secondary star (Verner et al. 2005; Duncan & White 2003; Soker 2007a). Duncan 
& White (2003) further suggested that most of the radio emission comes from a disk of a 
few arcsec, ~ 10 4 AU at the distance of rj Car, around the binary system. We accept the 
general picture suggested by Duncan & White (2003) and White et al. (2005), and try to 
explain the general properties of the radio emission. 

The 3.5 cm-radio emission after the 1998 minimum was lower than in the previous cycle 
(Duncan & White 2003; White et al. 2005). Duncan & White (2003) attributed this fading 
to a denser equatorial gas due to shell ejection at the 1998 periastron passage. The shell 
ejection idea goes back to Zanella et al. (1984). However, the shell ejection model fails to 
explain the properties of the X-ray minimum (Akashi et al. 2006; see discussion in Soker 
2007a), as well as the behavior of some visible emission lines (Nielsen et al. 2007). We will 
therefore examine other explanations, such as the secondary became cooler and changes in 
the magnetic properties of the primary wind. 

Our main goal is to understand the basic processes behind the variations in radio lu- 
minosity in order to learn on the interaction of the two winds, and to reveal signatures of 
accretion during the spectroscopic event. In addition, as will be evident (see also Soker 
2007a) the general flow structure and ionization and recombination rates we deal with are 
far too sensitive to the stellar and wind parameters to be modelled by a simple approach. 
Therefore, we will not try to fit each point in the radio and millimeter light curves, but 
only to approximately match the general behavior within a cycle, with the goal of finding 
the role of the primary wind and the role of the secondary ionizing radiation. We will also 
speculate on the reasons for the fast variations in the intensity of the millimeter emission 
near periastron (Abraham et al. 2005). 



2. SETTING THE FLOW STRUCTURE 

2.1. The Binary System 

The r] Car binary parameters used by us are as in the previous papers in this series 
(Soker 2005b; Akashi et al. 2006; Soker & Behar 2006; Soker 2007a), and are compiled by 
using results from several different papers (e.g., Ishibashi et al. 1999; Damineli et al. 2000; 
Corcoran et al. 2001, 2004; Hillier et al. 2001; Pittard & Corcoran 2002; Smith et al. 
2004; Verner et al. 2005). The assumed stellar masses are M l = 12OM , M 2 = 3OM , 
the eccentricity is e = 0.9, and orbital period 2024 days, hence the semi-major axis is 
a = 16.64 AU, and the orbital separation at periastron is r = 1.66 AU. The mass loss rates 
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are Mi = 3 x 10 4 M yr 1 and M 2 = 10 5 M yr 1 . The terminal wind speeds are taken to 
be i>i = 500 km s _1 and v 2 = 3000 km s _1 . 

The dependence of the primary wind properties on latitude (Smith et al. 2003) is not 
important here, because we are interested in the interaction winds region which occurs near 
the equatorial plane. High latitudes are less relevant to our goal. The properties of the 
primary wind segments that interact with the secondary wind are mainly constrained by 
the X-ray emission (Pittard & Corcoran 2002; Akashi et al. 2006), and are used here. The 
general flow structure is depicted in Figure 1. 

Verner et al. (2005) deduced the following secondary stellar properties: T 2 ~ 37, 200 K, 
L 2 ~ 9.3 x 10 5 L Q , R 2 ~ 23.6i? , v 2 = 2000 km s~\ and M 2 ~ 8.5 x 1O~ 6 M yr -1 . The 
properties of the secondary star used by us are discussed in previous papers of the series 
(Soker & Behar 2006; Soker 2007a), and are based on the results of Verner et al. (2005). We 
take R 2 ~ 2Oi? and T 2 ~ 36, 000 - 40, 000 K. 

Using the results of Schaerer & de Koter (1997) for ionizing flux from hot stars, we fit 
the rate of ionizing photons ihv > 13.6 eV) emitted by the secondary star per steradian by 

^ = N^ = L 2 ^ f 3.5 x 10 48 T 2 = 40, 000 K 
12 4tt 9 x 1O 5 L \ 2.3 x 10 48 T 2 = 36, 000 K 1 ' 

where a linear fit can be done for any effective temperature in the above range (see Soker 
2007a). 



2.2. Undisturbed Primary Wind 

We examine the fraction of the ionizing photons absorbed by the undisturbed (freely 
expanding) primary wind. Namely, taking the region from r 2 = r 2 o to r 2 = r 2 i to be the 
freely expanding primary wind, i.e., the secondary stellar wind and gravity do not deflect 
the primary wind in that region. Both r 20 and r 2i can depend on the angle a. Inside the 
low-density cone the density of the secondary wind is very low, and it can be neglected when 
calculating the recombination rate. The different variables that will be used are defined in 
Figures [1] and [2j 

The density of the primary wind as function of distance r 2 from the secondary is 

Mi 

Anv i(r 2 + r\ + 2r 2 r cos a) ' 

where r is the orbital separation and a the angle between the direction of r 2 (measured from 
the secondary) and the line joining the two stars. The total recombination rate per steradian 
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discontinuity 



Fig. 1. — Schematic drawing of the collision region of the two stellar winds and the definition 
of several quantities used in the paper. There is an axial symmetry around the line through 
the two stars. The two thick lines represent winds' stream lines. The two shock waves are 
drawn only in the lower half. The post-shock regions of the two winds are hatched. The 
region dominated by the low-density secondary wind, to the right of the contact discontinuity 
in the figure, is termed the low-density cone (or shock cone). The dashed line shows the 
accretion column which exists, according to the proposed model, only for ~ 70 — 80 days 
during the accretion period which corresponds to the X-ray minimum and the spectroscopic 
event. 
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along that direction is 



r-ii 



R(r 2 o,r 2i} a) = a B n e n p r 2 dr 2 = G(M 1 ,v 1 ) 



''20 



T2l 



''211 



( 2 , 2 , 1 ^dr 2 , (3) 

[r z + ri + 2r 2 r cos a)" 1 



where as is the recombination coefficient, fimn is the mean mass per particle in a fully 
ionized gas (we assume solar abundance), and 



G?(Mi,Ui) = 0.22a B 



Mi 



(4) 



The integral in equation ([3]) has an analytical solution that gives for the recombination 
rate per steradian 



R(r 20 ,r 2i ,a) 
G(M l}Vl ) 



2r 2 cos 2 a — r 2 + r cos a 



2 sin air 2 + r| + 2r 2 r cos a) 2r sin a 



tan" 



_i r 2 + r cos a 



r sin a 



''2i 



''20 



(5) 



Taking r 2 i —> oo and r 2 o = gives the absorption in the case where the primary wind is not 
disturbed at all 



R(0,oo,a) = G{M 1: v x ) 



Mi 



a cos a 

- 2r sin 3 a 2r sin 2 a 
Substituting typical values for G{M\,V\) in equation (jGj) we find 

2 

/ \ r \ 

R{0,oo,a) = 3.2xl0 47 



(6) 



-i 



a 



cos a 



-i 



3 x 1O- 4 M yr- 1 J VSOOkms- 1 ^ V5 AU J V s in 3 a sin 2 a 

(7) 

We note that the freely expanding primary wind absorbs radiation mainly in the direc- 
tion a > 140°, namely in directions toward the primary star. However, equation (jSJ) largely 
overestimates the absorption of the secondary ionizing radiation in those directions because 
the primary itself ionizes the material close to it. 



2.3. The Shocked Primary Wind 

The secondary wind 'cleans' the area behind the secondary (right side in Figured]) and 
compresses the primary wind along the contact discontinuity, increasing the recombination 
rate there. For the parameters used here the half opening angle of the wind-collision cone is 
cf) a ~ 60° (Akashi et al. 2006). This implies that even as the secondary approaches or recedes 
periastron the secondary fast wind will clean a large solid angle for ionization to propagate 



- 6- 




Fig. 2. — Schematic drawing of the collision region of the two stellar winds and the definition 
of several quantities used in the paper. The contact discontinuity is approximated as having 
hyperbolic shape. The wide arrows represent absorption of ionizing photons: The empty 
arrow is in the secondary wind region (low-density cone), where absorption is negligible. The 
short black arrow is in the primary wind postshock region, where most absorption occurs. 
The hatched arrow represent absorption in the freely expanding primary wind region. The 
figure is not to scale. For the parameters used here D g2 — 0.3r, and d p <C r 2s in the relevant 
regions. 
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Fig. 3. — The orbital plane and a cut through the primary wind shocked region; this region 
extends beyond the size of the figure, as drawn schematically at phase t = —0.06. The 
phases are marked at four points along the orbit. The elliptical orbit, the primary stellar 
size, the distance of the shocked primary wind region and its width, are drawn approximately 
to scale (but not the secondary star or the secondary wind shock region). The thick section 
of the orbit is approximately where accretion phase occurs according to the model. At 
phase t = 0.5 we mark the sources of the X-ray emission (the shocked secondary wind), of 
the ionizing radiation (the secondary photosphere). Two possible orbital orientations are 
shown. Position 1: The secondary is closer to us near periastron passage. Position-2: The 
primary is closer to us near periastron passage (Nielsen et al. 2007). In both cases the 
orbital plane is tilted by ~ 45° to the line of sight. 
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almost unattenuated at low latitudes, unless the fast wind is shut down for ~ 10 weeks as 
assumed here (Soker 2007a). The general structure along the orbit is shown at 4 epochs in 
Figure [3 

The fast orbital motion near periastron causes this low-density (shock) cone to wind 
(wrap) around the secondary, and therefore to increase absorption (Nielsen et al. 2007). This 
is important only very near to periastron passage, and it is much too complicated to deal 
with in this analytical paper. But most important, when this winding becomes important 
we expect an accretion phase to start (Soker 2005b). We will therefore, ignore winding of 
the wind interaction region around the binary system. 

The primary wind postshock region is unstable, and has a corrugated structure (Pittard 
& Corcoran 2002; Pittard et al. 1998). Magnetic fields within the primary wind will reduce 
the compression factor as we shall see below. Because of the unstable postshock flow structure 
and the unknown magnetic field intensity, our calculation here is a crude estimate, but still 
teaches us on the importance of the shock wave. The post-shocked primary wind flows in a 
thin shell along the contact discontinuity with a speed v^, as drawn schematically in Figures 
[Hand [2], where all variables used here are defined. Let the distance a ray from the secondary 
cuts in the thin shell be d p , such that the shell width (perpendicular to the flow) is d p cos5. 

The contact discontinuity shape is approximated by an hyperbola at distance D g 2 from 
the secondary at the stagnation point of the colliding winds, where the two winds momenta 
balance each other along the symmetry axis, and asymptotic angle <ft a = 60°. The velocity 
Vd of the post-shock primary wind parallel to the shock front is zero at the stagnation point, 
and increases to ~ v i at infinity. Mixing between the two winds, as a result of instabilities 
(the corrugated structure; Pittard & Corcoran 2002), will further accelerate the post-shock 
primary wind (Girard & Willson 1987). This happens because the post-shock secondary 
wind expands faster than the primary wind and because of its long cooling time it has a 
pressure gradient parallel to the contact discontinuity. Therefore, the value of Vd is hard to 
calculate in the analytical model. We therefore take this velocity to be 

Vd = v x q v cosil), (8) 

where q v is a parameter that is kept constant in the model. Increasing q v will have similar 
effect as reducing the density compression ratio (see below). We take q v = 0.5. 

The density of the postshock primary wind introduces large uncertainties, because it is 
strongly influenced by the magnetic field. The primary wind rapidly cools to a temperature 
of T p ~ 10 4 K, and it is compressed by the ram pressure of the slow wind to a density 
p p . Neglecting first the magnetic pressure in the post-shock region, p p is given by equating 
the thermal pressure kT p p p / prrin of the post-shock material with the ram pressure of the 
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primary wind pi(fwindi sin^) 2 , where ip is the angle between the slow wind speed and the 
shock front, and f w i nc n is the preshock speed of the primary wind relative to the stagnation 
point;. Because in most interesting orbital phases the orbital velocity is low, we take here 
^windi = Vi. For a post-shock temperature of T p ~ 10 4 K we find for the compression factor 
neglecting the magnetic field 

PA - 500 f ^— T YY— Y- 0) 



.Pi /thermal V500 km B"V V 0.5 

Such a high density contrast is seen in Figure 2 of Pittard & Corcoran (2002). However, 
because of magnetic fields we don't expect such a large compression behind the shock. 
Typical preshock magnetic pressure to ram pressure ratio can be 

rjB^Pm/pivf- 0.001 -0.1 (10) 

(e.g., Eichler & Usov 1993; Pittard & Dougherty 2006). The magnetic field component 
parallel to the shock is increased as it is compressed when the density is increased in the 
shock wave. For a random field we can take this component to contribute 1/3 to the preshock 
pressure. Equating the post-shock magnetic pressure to the wind's ram pressure we find the 
limit on the compression factor imposed by the magnetic field to be 

Wb \W Vo.002/ 1 0.5 ( ' 

The strong magnetic field can also smooth the strong corrugated structure seen in the sim- 
ulations presented by Pittard & Corcoran (2002). The presence of the magnetic field is 
another source of the large uncertainties involved in our calculation, and it can introduce 
large stochastic variations on short time scale and from cycle to cycle. We will use equation 
ffTTj) to calculate the density p p because (p p J P\)b < (pp/pi) thermal in the cases studied here. 
In Figure H] we plot p p /pi as function of the angle a. 

Mass conservation for the primary wind implies that the mass expelled within an angle 
(f) equals the flow in the postshock shell. This reads 

Mi(0) = -(1 — cos<p) = 27ir 2s sin aVdPp cos 5d p , (12) 

where r2 S is the corresponding distance of the contact discontinuity from the secondary. The 
quantity d p is calculated from equation (TT21 . The ratio <i p cos5/r 2s as function of a is plotted 
in Figure [51 

The recombination rate per steradian in the postshock shell is given by 

-Rshock = o> B (n e n p ) p rl s d p , (13) 
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Fig. 4. — The density compression factor in the shock of the primary wind for preshock 
magnetic pressure to ram pressure ratio of t)b = 0.002. In the entire range the density ratio 
is determined by equation (jTTi) . demonstrating the very important role of the (weak) magnetic 
field in the preshock primary wind. If equation ffTTl) gives p p / pi < 1 we set p p / pi = 1 
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Fig. 5. — The ratio of the postshock shell width d p cos 5 to r 2s for the standard parameters 
listed in section [3] (see Figure [2] for definitions of these quantities). 
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where, as before, subscript 'p' refers to postshock quantities (e.g., (n e n p ) p is the value of 
the electron number density times proton number density calculated for the postshock gas). 
Using equation (TJJ we can scale equation (TTHT) and write 



2.4. The Radio Emission 

The optical depth for radio emission (e.g., eq. 4.32 in Osterbrock 1989) is given by 

T \~ L35 / A \ 21 / n P \2 / 



T ~°' 1 \10*k) \cm) ll0 6 cm- 3 J \l00 AU ) ^ 

where h is the typical size of the region considered, and n e is the electron density. For the 
undisturbed primary wind and fully ionized gas 

n e = 7.56 x 10 5 ( ^ -) { ? \ " (-?*-) " cm- 3 . (16) 

^3 x 1O" 4 M yr- 1 / V500 km s- 1 / V 300 AU J v ; 

We see from equation ffl5l) that for A ~ 3 cm the optical depth is r > 1 at distances 
n < 300 AU. 

The contribution of the volume inner to the distance r\ ~ 300 AU, where r(3 cm) ~ 
1, is small because it is optically thick at these wavelengths. We therefore take V2i = 
300 AU, which is ~ 10 times larger than apastron distance, in equation (j3j) when calculating 
absorption of ionizing photons by the freely expanding primary wind. The lower limit in 
equation (jSJ) is r 2 o = r 2s . 

Most radio emission comes from much larger distances, up to ~ 10 4 AU (Duncan & 
White 2003; White et al. 2005). There we expect that the emission will come from dense 
clumps formed by instability of the outflowing primary wind. We therefore assume that 
many clumps within a distance of ~ 1" = 2300 AU, and density as the wind density is at 
r\ ~ 300 AU from the primary star, absorb a large fraction of the ionizing radiation escaping 
the binary system and they are responsible for a large fraction of the radio emission. They are 
optically thick at A ~ 3 cm. In addition, some ionizing radiation escape to larger distances, 
leading to no contribution to ionizing gas that emits in the radio band. We therefore use an 
efficiency parameter k r (X) to convert from maximum possible radio flux to observed radio 
flux. The maximum possible radio emission is for the case where all ionizing radiation ionizes 
optically thin gas to ~ 10 4 K in the observed region and in a steady state. k r increases with 
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decreasing wavelength with k r (3 cm) ~ 0.1. Higher intensities are expected for shorter 
wavelength, where optical depth is lower. Indeed, emission flux at A = 1.2 cm is ~ 2.5 times 
as bright as that for A = 3.5 cm in 2004 (White et al. 2005). 

We define the difference between ionizing photon rate emitted by the secondary, Ni2, 
and the total absorption 

N diS = N l2 - [R(r 2s , 300 AU, a) + R shock («)]. (17) 

This quantity can have negative values, which is meaningless. In calculating the emitted 
radio luminosity when ever N^g < we set it to zero. The radio emission at each orbital 
angle 9 is calculate according to the formula 

/>180° 

^radio = k r h2 , K I max(iV diff , 0) sin ada, (18) 

Ja=0 

where A is the conversion factor from ionizing photon rate (in equilibrium with recombina- 
tion) to radio power at 10 4 K (Osterbrock 1989). 



3. RESULTS 

3.1. The Basic Processes 

The parameters of the binary systems are given in section 2.1. The other parameters 
for our standard follows. The contact discontinuity hyperbola has a distance 

D g2 = 0.3r from the stagnation point and an asymptotic angle <p a = 60°. The velocity 
parameter in equation (JS]) is q v — 0.5, the magnetic pressure in the primary wind t]b = 0.002, 
the radio efficiency k r = 0.07, and the ionizing photon rate per steradian of the secondary 
star Ni2 = 2 x 10 48 s _1 sr _1 (eq. [1]). We could change the parameters in a way that will 
not change the results much. For example, we can increase Ni2 together with increasing a 
little the compression ratio in the shock (by lowering t]b) and reducing the radio efficiency. 
These parameters are sufficient for our goal of demonstrating the importance of absorption 
in the shocked primary wind, and connect this to the fast variation in millimeter emission 
near periastron (section H]). 

First we examine the importance of the absorption by the primary wind of the ionizing 
radiation emitted by the secondary star. In Figure we draw the absorption rate (assumed 
to be equal to the recombination rate) per steradian as function of the angle a as defined 
in Figure [2j a = 0° is away from the primary and a = 180° is toward the primary. The 
thick lines in the left panel show the absorption due to the shocked primary wind as given by 
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equation f|T4l) . and the thin lines show the absorption of the freely expanding primary wind 
as given by equation ([5]), from the shock position r 2 o = t 2s (q;) to a distance r 2 « = 300 AU. 
In the right panel we show the value of N^m (eq. [T7I) . the difference between ionizing photon 
rate emitted by the secondary, N diS (eq. [TTj), for N i2 = 2 x 10 48 s _1 sr _1 (dashed line in the 
figure), and the absorption rate. 

For a < 60° the secondary wind cleans a low-density cone, and no absorption occurs 
there; this region is not in the graph. (We don't consider the winding of this low-density 
cone around the primary near periastron passage, which will increase absorption in directions 
a < 60°; Nielsen et al. 2007.) The absorption by the free wind is overestimated for large 
values of a. The reason is that these directions pass close to the primary star, where the 
primary radiation ionizes a large volume around the primary. This is not considered in our 
calculation. 

The sharp change in the slope of the lines presenting the absorption by the shocked 
wind near a = 60° result from our constraint that the density compression ratio cannot be 
Pp/pi < 1, while close to the a = 60° equation (fTT|) gives a value < 1. The same change 
in behavior is imprinted in the plots of iVdiff in the right panel. The results close to the 
asymptotic lines of the hyperbolic shock wave are not trustable, but in any case this region 
plays a minor role. 

In Figure [7] we show the radio light curve for our standard case as calculated by equation 
ffT8"j) with k r = 0.07. We observe two orbital phases (and their two symmetric phases that 
are not shown) where a change in behavior occurs in the light curve. Near periastron, at 
phase< 0.011 (time< 22 day before or after periastron passage) corresponding to orbital 
angle 9 < 95°, the radio flux is constant. This is the region where all ionizing photons in 
our model are absorbed by the primary wind, beside those emitted into the low-density cone 
a < 60°. This region is not modelled correctly by us for two reasons. First, the low-density 
cone is winding around the binary system (Nielsen et al. 2007) as its orbital velocity is large 
(~ 300 — 400 km s -1 ) in this part of the orbit. Second, and more important, at this phase 
an accretion process is likely to start, where the secondary accretes from the primary wind 
(see section H]) and the low-density cone does not exist for ~ 10 weeks. 

The second change in behavior occurs at phase~ 0.094 (~ 191 day) corresponding to 
orbital angle ~ 152°. As can be seen from the right panel of Figure El this occurs when the 
absorption in the directions 80° < a < 150° drops below the secondary ionizing photon flux. 
For other parameters the change will occur at a different orbital phase. 
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Fig. 6. — Absorption rate and emission rate per steradian, for the standard parameters 
given at the beginning of section 3, given every 20 degrees of the orbital angle from 6 = 
(periastron) to 9 = 180° (apastron). Left panel: Absorption rate per steradian of ionizing 
photons by the primary wind, assumed to be equal to the recombination rate. The absorption 
rate of the shocked primary wind as given by equation (TH1) are drawn by thick lines. The 
absorption rate of the undisturbed primary wind, as given by equation ([5]), from the shock 
position r2o = r2 S (u) to a distance r2% = 300 AU are drawn by the thin lines. Right panel: 
The difference between the secondary ionizing rate -/Vj 2 = 2 x 10 48 s _1 sr _1 (dashed line) 
and the total absorption, as given by equation f[T7|) . t stands for the phase, with t = at 
periastron and t = 0.5 at apastron. 
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Fig. 7. — The calculated radio-flux light curve. F = L radio /47r.D 2 , where L Tadio is given by 
equation (fT8l) (with k r = 0.07), and D = 2.3 kpc is the distance to T) Car. The lower line is 
for the standard parameters (given at the beginning of section 3). The two other lines are for 
higher secondary ionizing photon rate, as marked near each line in units of ionizing photons 
per second per steradian. t stands for the phase, with t = at periastron and t = 0.5 at 
apastron. 
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3.2. Periastron Passage 



According to the accretion model the deep X-ray minimum is assumed to result from 
the collapse of the collision region of the two winds onto the secondary star (Soker 2005b; 
Akashi et al. 2006). This process is assumed to start ~ — 20 day before periastron passage, 
to last ~ 10 weeks, and shut down the secondary wind, hence the main X-ray source. During 
that period the secondary accretes mass from the primary wind at a rate of ~ 10~ 6 M© yr _1 . 
Akashi et al. (2006) showed that this assumption provides a phenomenological description 
of the X-ray behavior around the minimum. We examine the implications of the accretion 
process to the model here and in section 4. 

After the secondary wind is shut down there will be no low-density cone behind the 
secondary any more. The accretion process results in a high density region behind the 
secondary star, the accretion column (dashed line in Figured]). No radiation will escape in 
that direction. However, there is no shocked primary wind any more, and away from the 
equatorial plane the absorption will be that of the freely expanding primary wind according 
to equation (J7|). For our standard case we find that perpendicular to the line joining the 
two stars, a = 90°, and at apastron r = 1.66 AU, the absorption rate is -R(0,oo,90°) = 
1.5 x 10 48 s _1 sr _1 . Since R(0, oo, 90°) < N i2 , we expect ionizing radiation to escape to large 
distance. Therefore, the ionizing radiation is not expected to drop to zero even during the 
~ 70 day long accretion phase. 

But even if ionizing radiation drops to zero, the radio emission will not follow. The 
reason is that the recombination time (for our standard parameters) 



at relevant distances is longer than the accretion phase. Many dense regions which give more 
radio emission will recombine on a shorter time scale. This holds along the entire orbit. The 
recombination time at a distance of ~ 1" ~ 2300 AU is longer than the orbital period, and 
must be considered in trying to exactly fit the light curve. But this effect is much more 
important near periastron passages, when changes are fast and the wind segments as close 
as few hundred AU from the system don't have time to recombine. 

To summarize this section, it is possible that when the system is near periastron there 
is still some ionizing flux reaching large distances. But even if not, the accretion phase last 
only ~ 10 weeks, and the relevant regions of the wind have no time to fully recombine. 




(19) 
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3.3. The Role of the Different Variables 

Because the ionizing radiation that escape to large distances is the difference between 
the emitted and absorbed ionizing radiation, both the radio emission (section 3.1) and the 
highly ionized emission lines in the UV and visible bands (Soker 2007a) are very sensitive to 
the parameters of the system. 

The ionizing radiation rate Ni2 has an obvious influence. In Figure [7] we present the 
results for increasing A^2 by a factor of 1.15 to Ni2 = 2.3 x 10 48 s _1 sr _1 , and by a factor of 
1.5 to N i2 = 3 x 10 48 s -1 sr -1 . We see how the resulting radio emission at maximum increases 
by a larger factor than the increase in N i2 . Near periastron the radio emission results from 
region ionizes by the ionizing photons escaping only from the low-density cone, and hence it 
increases by the same factor as A^2 was increased. As discussed above, for ~ 10 week very 
close to periastron passage an accretion process is likely to occur, and the low-density cone 
does not exist. 

The role of the primary mass loss rate and velocity is clearly evident from the role of 
G(Mi,vi) ex M\jv\ in equations ([3])- (Q and (TH|) . The absorption rate is sensitive to these 
parameters, and therefore the ionizing radiation escaping to large distances. Changing these 
parameters not only affect the absorption by the wind, but also change the opening angle 
of the low-density cone. For the parameters used here the dependance is approximately 
(Eichler & Usov, 1993) <p a oc (M 2 f 2 /MiWi) 1/3 . The radiation escaping through the cone is a 
fraction (1 — cos0 a )/2 of the total secondary ionizing radiation. The sense of the changes in 
the radio luminosity is opposite to that of the primary mass loss rate, i.e., a higher value of 
Mi results in a lower radio luminosity. 

The primary wind speed, on the other hand, might have an opposite role near apastron 
and periastron. Say the primary wind decreases by a factor of ~ 1.5 to v\ ~ 330 km s~ 4 . 
This implies a higher density in the wind. For our standard case we find that for the entire 
orbit the total absorption is above the emission rate, such that the only ionizing radiation 
escape from the low-density cone. Using the formula for the low-density cone opening angle 
as given by Eichler & Usov (1993) we find that it increases from <ft a = 61° to <p a = 68°. The 
fraction of secondary ionizing photons that escape to large distances increases from ~ 0.26 
to 0.32. Remembering that near periastron most of the ionizing radiation escapes from the 
cone (but not during the ~ 10 weeks of the accretion phase very close to periastron passage, 
when the low-density cone does not exist), we expect an increase of radio emission by a factor 
~ 20% near periastron passage, but a decrease in radio emission near apastron, compared 
to our standard case. 

Another important parameter is the density compression ratio in the shock. Reducing 
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the magnetic pressure tjb in the preshock wind will allow larger compression, i.e., larger 
Pp/pi in equation ffTTj) . This results in higher absorption and lower observed radio flux. In 
Figure [H] we show our standard case plotted for the 1992-8 cycle, and a case of r] B = 0.001, 
with all other parameters unchanged, plotted on the 1998-2003 cycle. The observed radio 
light curve is from White et al. (2005). We don't try to fully fit the radio observations, but 
just to show the effect of varying the value rjs, and by that emphasizing its importance. 

To further illustrate the effect of some variables, in Figure M we repeat the presentation 
of Figure E] but we changed the ionizing flux from N i2 = 2 x 10 48 s _1 sr _1 to N i2 = 3.5 x 
10 48 s _1 sr _1 , and the efficiency parameter from k r = 0.07 to k r = 0.04. The values of tjb are 
given on the lines. 

4. FAST VARIATIONS NEAR PERIASTRON PASSAGE 

The accretion process near periastron passage might lead to the formation of an accretion 
disk for a very short time (Soker 2003, 2005a), and might even lead to a transient launching 
of two opposite jets, as suggested theoretically (Soker 2005a; Akashi et al. 2006), and might 
have been observed (Behar et al. 2006). Basically, in a steady state Bondi-Hoyle type 
accretion flow the accreted mass has not enough specific angular momentum to form an 
accretion disk. However, stochastically accreted blobs at the onset of the accretion phase 
might lead to the formation of a transient accretion disk, and possibly to two jets. 

Because the 3.5 cm observations (White et al. 2005) don't sample the time near peri- 
astron, we consider the 7 mm and 1.3 mm light curves reported by Abraham et al. (2005). 
In Figure [TD] we compare the temporal evolution of the X-ray flux (Corcoran 2005) with the 
7 mm flux and the 1.3 mm flux from Abraham et al. (2005) near periastron. Zero time is 
taken at June 29, 2003. The scales and vertical displacements are not the same for the dif- 
ferent graphs, as we are only interested in the temporal behavior. In Figure [TOj the 1.3 mm 
flux varies between 5 Jy and 23 Jy, and the 7 mm flux varies between 0.7 Jy and 3.4 Jy. The 
1.3 mm flux is much larger that the 7 mm flux as a result of contribution from hot dust to 
the 1.3 mm flux (Cox et al. 1995), and because optical depth is lower at 1.3 mm. 

There is some correlation between the X-ray and millimeter emission behavior. This can 
be understood as follows. From Figure [6] we learn that near periastron all ionizing radiation 
escape from the opening formed by the secondary wind, i.e., with in angle a < 60° (but 
not during the accretion phase). We ignored in our treatment the fast orbital motion near 
periastron that causes this opening cone to wind around the secondary, therefore increasing 
absorption (Nielsen et al. 2007). This winding of the winds-interaction region is not enough 
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Fig. 8. — The calculated radio- flux light curve and the observed one from White et al. (2005). 
F = L ra dio/47r-D 2 , where L ra dio is given by equation ( fT8l) (with k r = 0.07), and D = 2.3 kpc 
is the distance to rj Car. In the 1992-1998 period we use our standard parameters, while in 
the 1998-2003 cycle we take rj B = 0.001 instead of r} B = 0.002. We don't try to fully fit the 
radio light curve, but just to demonstrate that the parameters can fit the observation, and 
emphasize the role of the weak magnetic field in the primary wind. 
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Fig. 9. — Like FigureEl but the ionizing flux is increased by a factor of 1.75, and the efficiency 
of producing radio emission reduced from k r = 0.07 to k r = 0.04. The value of r\ B is given 
for each curve. 
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Fig. 10. — X-ray (Corcoran 2005) and millimeter (Abraham et al. 2005) emission from 77 Car 
at the beginning of the spectroscopic event, drawn to emphasize the peaks and deeps in the 
emission curve. The three graphs don't have the same flux scale and are shifted vertically 
arbitrarily. Zero time is at June 29, 2003. 
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to account for the X-ray decline, and for that the accretion model was proposed (Soker 
2005b; Akashi et al. 2006). The accretion process is likely to be stochastic, especially at 
the beginning. At the beginning we expect blobs of material to be accreted (Soker 2005b). 
This accretion is assumed to shut down the secondary wind. The stochastic accretion at the 
beginning of the accretion phase is likely to lead to a few weeks of sporadic ejection of the 
secondary wind, mainly perpendicular to the orbital plane. The sporadic fast ejecta from 
the secondary is the explanation of the high X-ray peaks observed just before the long X-ray 
minimum (Soker 2005b; Akashi et al. 2006). 

The sporadic ejection of the secondary wind will increase the radio emission as well. 
When dense blobs are accreted they will absorb more of the secondary radiation simultane- 
ously with the reduction in the secondary wind and X-ray emission. When the secondary 
wind is intensified and lead to X-ray flare, it also cleans a region around the primary, allowing 
more ionizing radiation to escape. This will increase the radio emission. 

The accreted matter might also form a blanket at a distance of a few solar radii, or 
~ 0.1 — 0.3i? 2 , above the secondary (Soker 2007a) . This reduces the effective temperature 
and as a result of that the ionizing flux. As evidence from Figure El the ionizing rate is close 
to the absorption rate. This implies that any small change in the ionizing flux can lead to 
a large change in the difference between emission and absorption. The last peak in 7 mm 
emission (a small knee is seen also in the 1.3 mm light curve at that time) that corresponds 
to no X-ray flare, might be a short heating episode of the secondary surface. 

X-ray flares exist along the entire orbit, but they are much stronger near periastron 
passages (Corcoran 2005). The X-ray variations along most of the orbit result from stochastic 
variations in the winds properties and instabilities in the winds collision process. According 
to the accretion model (Soker 2005b; Akashi et al. 2006) the larger variations near periastron 
passages result from instabilities in the accretion process. We propose the same type of 
behavior for the radio emission. Namely, large variations in the free-free radio emission 
near periastron passages result from instabilities in the accretion process that cause large 
variations in the escaping ionizing radiation. 

5. SUMMARY 

Our main goal was to explore the effect of the interaction of the two winds in rj Car binary 
system on the radio emission. The general flow structure and ionization and recombination 
rates we deal with are far too sensitive to the different parameters for us to build an exact 
model. Our results use, and support, previous claims that the 3.5 cm radio emission from 
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rj Carinae can be explained by free-free (bremsstrahlung) emission from ionized gas (Cox et 
al. 1995; Duncan & White 2003), and the ionization source is mainly the binary secondary 
star (Verner et al. 2005; Duncan & White 2003; Soker 2007a). What we have shown here is 
that the absorption of ionizing photons by the shocked primary wind must be considered in 
addition to the absorption by the freely flowing wind (Soker 2007a). 

The radio emission is likely to come from dense regions at distances of hundreds to 
thousands AU from the binary system. Therefore, the observed radio flux depends sensi- 
tively on the unknown distribution of mass in these region. We parameterized the radio 
emission efficiency of the ionized gas by a parameter k r . We found that k r ~ 0.07 in our 
phenomenological model. 

The optical depth of the primary wind to 3 cm radio emission at distances smaller than 
~ 300 AU from the system is large (eq. [15]) . and most radio emission comes from larger 
distances. The optical depth is smaller at smaller wavelengths. Therefore, the exact radio 
light curve will be different at different wavelengths. This different behavior at different 
wavelengths (2, 6, 20 cm) is seen in the very eccentric massive binary star system WR 140 
(HD193793; White & Becker 1995). However, the mass loss rates in WR 140 are much lower 
than those in rj Car, and no direct comparison can be made. 

The ionizing photon rate reaching large distances from the secondary is the difference 
between the emitted rate by the secondary star (Verner et al. 2005; Duncan & White 
2003; Soker 2007a), with some contribution from the primary, and the absorption rate. The 
influence of the ionizing photon rate is explored in Figure [71 

The absorption rate is sensitive to the primary wind mass loss rate and velocity, and to 
the compression of gas in the shock. As the primary stellar wind collides with the secondary 
wind it passes through a shock wave, cools to low temperatures (~ 10 4 K), and is compressed 
by the ram pressure of the winds until its magnetic pressure equals the ram pressure. Our 
main new finding is that the compressed post-shock primary wind has a crucial influence 
on the absorption rate, hence on the ionization state of the circumbinary material, and 
consequently on the observed radio flux. The absorption rate by the freely-flowing primary 
wind and by the shocked primary wind is presented in Figure The magnetic pressure 
determines the compression factor of the postshock primary wind. This is parameterized by 
t)b (eq. [TO]) , which is the ratio of the preshock magnetic pressure to the ram pressure of the 
primary wind. The influence of T]b on the observed radio flux is demonstrated in Figure [HJ 

The role of the magnetic fields in the primary wind might be far beyond determining 
the compression factor. Soker (2007b) speculated that the 20 years long Great Eruption of 
rj Car that occurred 160 years ago was triggered by magnetic activity in the inner convective 
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region of 77 Car. This activity expelled the outer radiative zone which is now the bipolar 
nebulae, the Homunculus. 

The winds properties determine also the opening angle of the low-density (shock) cone, 
(p a (Figute|2J). Close to periastron passage most, or even all, of the ionizing photons escape 
through the low-density cone (again, for ~ 10 weeks very close to periastron passage the 
low-density cone does not exist in our model because of the accretion process). Therefore, 
the opening angle of the cone is another source of sensitivity of the radio emission on the 
winds properties. 

In section 4 we commented on the fast variation in the 7 mm and 1.3 mm emission 
reported by Abraham et al. (2005). We accept the proposed accretion model for the spec- 
troscopic event (Soker 2005b; Akashi et al. 2006; Soker & Behar 2006; Soker 2007a). In 
that model for ~ 10 week near periastron passage the secondary star accretes mass from the 
primary wind instead of blowing its own wind. The accretion phase implies the following. 
(1) There is no low-density cone any more. (2) A high density accretion column (see fig- 
ured]) absorbs all radiation to the direction opposite that of the primary (a ~ 0); (3) The 
accreted mass can form a blanket around the secondary, reducing its effective temperature 
and the ionizing flux from its surface (Soker 2007a). (4) At the onset of the accretion phase 
the accretion process can be stochastic and the accretion material posses a relatively large 
specific angular momentum. This process can lead to the sporadic launching of a collimated 
secondary wind. 

The sporadic launching of the secondary wind can lead to X-ray flares (peaks in X-ray 
emission), and reduces absorption of ionizing radiation by clearing a path through the dense 
primary wind. This is our explanation to a possible correlation between X-ray peaks and 
millimeter emission peaks (figure [TUT) . The next spectroscopic event is predicted to occur in 
January 2009. We predict that near that time, i.e., fall of 2008 to January 2009, there will be 
some correlations between the X-ray and radio intensity. We also predict that signatures of 
collimated outflow might be detected during these 'flares'. For example, strong blue-shifted 
X-ray lines (Behar et al. 2006). 

Some processes were not included in our calculation of the observed radio flux. Future 
more sophisticated models will have to include these effects. In particular we did not consider 
the following processes. 

1. We did not consider the winding (wrapping) of the low-density cone around the binary 
system near periastron passage (Nielsen et al. 2007). This process affects results near 
periastron passage, but not too close, when our model assumes that the low-density 
cone does not exist because of the accretion phase. For most of the rest of the time 
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the winding has a small effect. 

2. The contribution of the primary star to the ionizing flux was not considered here. 

3. The dependence of the primary wind properties on latitude (Smith et al. 2003) was 
ignored in our treatment. 

4. The long recombination time at large distances (eq. [T9~j) . implies that in a more so- 
phisticated model the response time of the ionized gas to the ionizing radiation should 
also be considered. 

We thank Diego Falceta Goncalves and an anonymous referee for useful comments. 
This research was supported by a grant from the Asher Space Research Institute at the 
Technion. 
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Table 1: Binary Parameters 



Symbol 


Meaning 


Typical value 


r 


Binary separation 


1.66-31.6 AU 


e 


Eccentricity 


0.9 


Mi 


Mass of primary 


12OM 


M 2 


Mass of secondary 


3OM 


Mi 


Mass loss rate by primary 


3 x 1O" 4 M yr- 1 


M 2 


Mass loss rate by secondary 


1(T 5 M yr" 1 


«1 


Velocity of primary wind (relative to primary) 


500 km s" 1 


^2 


Velocity of secondary wind (relative to secondary) 


3000 km s- 1 


^windl 


Velocity of primary wind relative to secondary 


- 500 km s" 1 


i?2 


Stellar radius of secondary 


~ 0.1 AU 


^2 


Secondary Luminosity 


9 x 1O 5 L 


T 2 


Secondary effective temperature 


36,000 -40,000 K 
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Table 2: Geometrical and Other Parameters 



Symbol 


Meaning 


Typical value 


Defined 


Di 


Distance of stagnation point from primary D\ = r — D g2 


0.7r 


Fig. 1 


D 9 2 


Distance of stagnation point from secondary 


0.3r 


Fig. 1 




Asymptotic angle of contact discontinuity 


60° 


Fig. 1 


<P 


Angle from primary to a point on contact discontinuity 


0-60° 


Fig. 2 


a 


Angle from secondary to a point on contact discontinuity 


- 180° 


Fig. 2 


i> 


Angle between primary wind and contact discontinuity 


0-90° 


Fig. 2 


5 


Angle between secondary wind and contact discontinuity 


0-90° 


Fig. 2 





Orbital angle (6 = at periastron) 


-180° - 180° 




t 


Orbital phase 


- 1 (0 - 2024 day) 






Distance from primary to contact discontinuity 


~ 1— hundreds AU 


Fig. 2 


T 2s 


Distance from secondary to contact discontinuity 


~ 1— hundreds AU 


Fig. 2 


r 20 


Inner boundary of recombining undisturbed primary wind 


> 300 AU 


eq. © 




Outer boundary of recombining undisturbed primary wind 


~ r 2 s 


eq. © 


r 2 


Distance from secondary 


0.5— hundreds AU 




d p 


Length a ray from secondary cut in the thin shell 


~ 0.2- hundreds AU 


Fig. 2 


Vd 


Post-shocked primary wind speed in the thin shell 


- 250 km s" 1 


Fig. 2 


Qv 


Parameter determining the value of Vd by eq. (jHJ) 


0.5 


eq. © 


Pi 


Primary wind density 


Mi/4-KVxrl 


eq. © 


P P 


Density of post shock primary wind in the conical shell 


(1 - 40) x pi 


eq. (HU 


Vb 


Preshock magnetic pressure to ram pressure ratio 


0.001-0.1 


eq. (HH 


R 


Recombination rate in the undisturbed primary wind 


- 10 48 - lO^sr" 1 s" 1 


eq. © 


^shock 


Recombination rate in the conical shell 


- 10 47 - 2 x lO^sr" 1 s" 1 


eq. (0© 


n p 


Proton density 


10 5 - 10 12 cm~ 3 




n e 


electron density 


10 5 - 10 12 cm" 3 




N i2 


Rate of emission of ionizing photons by the secondary 


2 - 3.5 x 10 48 s" 1 


eq. © 




Efficiency of radio emission 


0.01-0.1 


eq. M 



